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Abstract. We calculate the change in the superconducting transition temperature 
T c of MgP>2 caused by interband nonmagnetic impurity scattering using the Eliashberg 
theory for the two-band model of this compound. Much slower rate of T c suppression 
is obtained compared to the prediction based on the BCS treatment of the two-band 
model which ignores renormalization and damping associated with the electron-phonon 
interaction. Hence, the interband impurity scattering rates deduced from experiments 
on MgB2 using the formula which results from the BCS approach to the two-band 
model are underestimated. We generalize the BCS treatment of the two-band model to 
include renormalization effects of the electron-phonon interaction and find an excellent 
agreement with the full strong coupling calculation. 
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1. Introduction 

There is a large body of experimental [1-14] and theoretical (for a review see [15]) 
evidence that MgB 2 is a multiband superconductor which is well described by an 
effective two-band model [16]. In the case of a multiband superconductor one expects 
that the superconducting transition temperature T c is reduced by the interband 
nonmagnetic (i.e. normal) impurity scattering in analogy to the effect of such scattering 
on anisotropic single band superconductors [17,18]. Several years before the discovery 
of superconductivity in MgB 2 the problem of impurity scattering in a multiband 
superconductor was examined in detail by Golubov and Mazin [19] using the weak 
coupling BCS-type treatment of the pairing interaction. They obtained an equation for 
the change in T c with the interband impurity scattering rate which is analogous to the 
Abrikosov-Gor'kov formula for the T c - suppression by paramagnetic impurity scattering 
in ordinary superconductors. The BCS-type treatment of Ref. [19] predicts that the 
T c is reduced by about 40 % for the interband scattering rate comparable to ksT c . 
For MgB 2 that would imply a drop in T c from 39 K to about 25 K for the interband 
impurity scattering rate Y = l/(2r) = 7/2 of about 1.7 meV. Thus, it was thought that 
observation of T c suppression with increasing disorder would provide the final evidence 
for the two-band model of MgB 2 . 

Experimentally, however, the situation appears to be more complicated. On one 
hand, as pointed out in [20] , the transition temperatures of different samples of MgB 2 
are rather insensitive to their respective residual resistivities : the T c s of samples with 
residual resistivities in the range from 0.4 to 30 /zffcm differ by at most 5%. On the other 
hand, irradiation of a polycrystalline sample of MgB 2 by fast neutrons led to an increase 
of residual resistivity and reduction in T c by as much as 20% [13]. The apparent lack 
of correlation between T c and residual resistivity in unirradiated samples was explained 
[20] by very small values of the interband impurity scattering matrix elements because of 
the particular electronic structure of MgB 2 so that the DC transport in this compound 
at low temperatures is primarily determined by intraband scattering which does not 
affect T c [19], while very weak interband scattering leads to no significant change in T c . 
The arguments in [20] apply to common substitutional impurities in MgB 2 which do not 
distort the lattice and subsequently Erwin and Mazin [21] proposed that substituting 
Mg with Al and/or Na would produce lattice distortions that could lead to large enough 
interband impurity scattering rates to cause a reduction of T c by a couple of degrees as 
predicted theoretically in [19]. Presumably the irradiation by fast neutrons generates 
enough lattice distortions to cause a 20% drop in T c [13]. 

Nevertheless, the break junction tunneling experiments on MgB 2 [4] clearly indicate 
that the interband impurity scattering is significant even in undoped and unirradiated 
samples. Namely, the only justification for using the equations of McMillan tunneling 
model for proximity effect [22] in analyzing the break junction data on MgB 2 is provided 
by the work of Schopohl and Scharnberg [23] on tunneling density of states of a 
disordered two-band superconductor. The fact that in the latter case the equations 
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have the form identical to those of the McMillan tunneling model for proximity effect 
is a pure accident as is evident from the entirely different meaning of the quasiparticle 
scattering rates in the two cases. The interband scattering rates used to fit the tunneling 
data [4] were at least as large as those predicted for Al/Na doped MgB 2 (Ts were 
in the range from 1 to 4 meV), but the T c of the material was reported to be 39 
K - close to the maximum value for MgB 2 of 39.4 K. A possible solution to this 
contradiction is that the weak coupling BCS-type treatment of impurity scattering 
in a multiband superconductor used in [19] is not quantitatively accurate for MgB 2 . 
The calculated electron-phonon interactions in MgB 2 [16] indicate that it is a medium- 
to-strong coupling superconductor (the largest calculated electron-phonon parameter 
X aa for o"-band electrons is comparable to the one in Nb) and renormalization and 
damping effects could play an important role in determining the rate of T c suppression 
by interband impurity scattering. 

In section 2 we solve the Eliashberg equations for a two-band superconductor with 
nonmagnetic impurity scattering and calculate the transition temperature as a function 
of the impurity scattering rate using realistic interaction parameters for MgB 2 [16,24]. 
We find that the T c is suppressed by interband scattering at much slower rate than what 
was obtained using the BCS treatment in Ref. [19]. In the same section we present the 
functional derivatives ST c /Sa 2 Fij, i,j — a, 7r [25] for several representative impurity 
interband scattering rates which show how the sensitivity of T c to various electron- 
phonon couplings changes with impurity scattering. In section 3 we generalize the BCS 
approach to include the renormalization caused by the electron-phonon interaction by 
extending the well known 9-9 model [26] to the two-band case. The numerical solution 
of such a model is found to be in excellent agreement with the full strong coupling 
calculation. In section 4 we give a summary. 



2. Strong coupling calculation 

2.1. Formalism 

The Eliashberg equations for T c of a superconductor with several isotropic bands 
i = 1,2,... which include nonmagnetic impurity scattering described by the Born 
approximation are [25,26] 

0*(n) = 0^)+E^^TV (1) 



+ °° 6-(m) 
<P» = nT c Y: E Mn - m) - ^0(E F - \ Um \)] (2) 

j m=-oo \UJ m \/jj{m) 



u n Zi(n) = uj n Z°(n) + E W—y- 1 ] » ( 3 ) 

+ °° CO 
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Here (f>i(n) is the pairing self-energy in the band i at the Matsubara frequency u n = 
7rT c (2n— 1) and Z^n) is the corresponding renormalization function (for a review of the 
Eliashberg theory of superconducting T c see [26]). The part <p®(n), Eq. (2), results from 
the intrabend and interband electron-phonon and screened Coulomb interactions, while 
the second term in Eq. (1) represents the impurity scattering contribution. In the same 
way, Zf(n) is the contribution to the renormalization function from the intraband and 
interband electron-phonon interaction and the second term in Eq. (3) gives the impurity 
contribution to Zi(n). The cutoff Ep on the sums over the Matsubara frequencies u m in 
Eq. (2) is initially taken to be large enough so that the Coulomb repulsion parameters 
are given by /i^ = VfiNj, where is the Fermi surface averaged screened Coulomb 
matrix element between the states in the bands i and j (Vfi = V£) and Nj is the 
Fermi surface density of states in band j [25] . The electron-phonon coupling functions 
a 2 Fij(Q) = a 2 fij(fl)Nj (a 2 fij(fl) = a 2 fji(Q)) enter via parameters A 



u 



n — m) 



r+oo 2Q 

Xijin -m)= dQa F ij (ti) —— J - , 2 , (5) 

and the impurity scattering rates ^ = l/r^ are given by (we use the units in which 
h—1 and the Boltzmann's constant ks—l) 

-!- = UimpTcNj | Vij | 2 (6) 

ZTij 

where n imp is the concentration of nonmagnetic impurities and is the Fermi surface 
averaged matrix element of the change in the lattice potential caused by an impurity 
between the states in the bands i and j. Clearly, 7^/7^ = Nj/N k = \j/\k, where 
Xij = Ay(0) (see equation (5)). 

In principle, Eqs. (l)-(4) have the form of an eigenvalue problem of a temperature 
dependent matrix with eigenvector 0, and T c is determined as the heighest temperature 
at which the largest eigenvalue of the matrix is one. However, before such solution is 
attempted one can simplify the problem further. First, by introducing the gap function 
as renormalized pairing self-energy Aj(n) = <f>i(n) j 'Zi(n) one can eliminate the intraband 
impurity scattering from the problem by combining Eqs. (1) and (3) 

Mn) U{n) + £ -L*^ = m + ]T-L^M (7) 

\ j^i ZT ij \ u n\ J j^i zr ij l^nl 

where 0°(n) is given by Eq. (2) with (ftj(m)/Zj(m) replaced by Aj(m). Next, the 
eigenvalue problem can be symmetrized by defining 

-Ai(n) 



Ui{n) = JNi- 



\ 



MZ?(n)+-£±- (8) 



together with Af 3 -(n — m) — ^ r N~jW j \ ij (n - m), n\- = ^Ni/Nj^ and 1/(2t£) = 
^Ni/Nj/(2Tij). With these definitions Eq. (7) reduces to 

00 XfAn — m) — n'iMEp - \u rn \) + ~A-r(l - 5 i3 )8 nm 
u l {n) = e(T)Y J E * T 7 , = = 7j j = f «*("») , (9) 

j rn=-oo ^UJ^Z^^UJ^Z^m) 
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with 




(10) 



and e(T)=l when T = T c . Finally, the size of the matrix which has do be diagonalized 
can be reduced by cutting off the Matsubara sums in (9) at a smaller energy oo c which 
is still large enough so that Z'An) ~1 for \uj n \ > u c ; hence, u c has to be at least 5-10 
times the maximum phonon energy Q m in various spectral functions a 2 Fij(Q) and much 
larger than the largest band off-diagonal l/2r^. The reduction in cutoff from Ep to 
oj c is accompanied by replacement of the Coulomb repulsion parameters //£■ in (9) with 
fj,*Au c ) where the matrix (in band indices) p,*(iv c ) is related to matrix fi s by [25] 



2.2. Numerical Results 

We solved Eqs. (9)-(ll) using the spectral functions a 2 F aa , a 2 F an , a 2 F nn and ct 2 F na for 
MgB2 obtained from the first principle electronic structure calculations and presented 
in [16]. The corresponding coupling parameters given by Eq. (5) with u n — uj m = are 
\ acr = 1.017, \ an = 0.212, \ nn = 0.446 and K a = 0.155. Since a 2 i^(ft) = a 2 /ij(^)-Nj, 
with a 2 fij(Vi) = a 2 fji(il), these values of A-parameters fix the ratio of the partial band 
densities of states N^/N^ = X a7r /X n(7 at 1.37. That fixes the ratio 7 <T7r /7 7rcr (see Eq. (6)) 
and we chose 7 7r<T as the independent scattering parameter. 

To minimize the effect of changes in the number N c = [uo c / (27rT c )+0.5] of Matsubara 
frequencies ([•••] denotes the integer part) on our numerical results as T c is reduced by 
increased interband impurity scattering rate we had to take the cutoff u c to be at least 
10 times the maximum phonon energy Q m . With u c fixed at 1000 meV the Coulomb 
repulsion parameters fi^ n , fj,* n , p0^ c were determined as follows. Choi et al. [24] 
calculated the ratios of the screened Coulomb repulsion parameters for MgB 2 to be \i aa 
: fi nn : fi U7T : fi na = 1.75 : 2.04 : 1.61 : 1.00. Since /i^ = V£Nj, the ratio Ha-wln-wo- 
= 1.61 implies that in their calculation N^/N^ = 1.61, which is considerably higher 
than the value of 1.37 found in [16] and adopted by us in this work by our choice of 
the electron-phonon coupling spectra. Leaving aside the reasons for such a discrepancy 
in N n /N a between the two sets of electronic structure calculations we use the ratios of 
the /x-values calculated in [24] to extract from them the ratios of the screened Coulomb 
matrix elements: V£ a : V£ a = 1.75 : 1.00, V£ n : V an = 2.04 : 1.61 and, because V% a = 
V^, V£ a : V£ n = 1.75 : 1.267. These could be combined with N n /N a = 1.37 to produce 
the ratios fiaa/Hmr = 1.01, Haaj H-an = 1-28 and fiaa/Hwa = 1-75 leaving the single fitting 
parameter \i aa once a choice is made for the initial cutoff Ep (see Eq. (11)). Since Ep 
is on the order of the largest electronic energy scale in the problem, we took E F to be 
equal to the 7r-bandwidth of 15 eV [27] and fitted \i aa in the four equations implied by 
the 2x2 matrix equation (11) 




(11) 



/C(^c) 



t^a-K f^ira ) hi 



Ep 



/D, 



(12) 
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E F - 



l^mr ~l" {f^a-af^wn A* unfair a) hi ' 



, (14) 
(15) 

PjJT / E J? \ ^ 

-0=1 + (/W + A^Tr) ^ (fiaaUnn ~ Hair Una) hi , (16) 

to the experimental T c o of 39.4 K for the case of no impurity scattering. The results were 
AW = 0.848234 with aCK) = 0.225995, aOc) = 0.225010 and aCK) = aOc) = 
0.067148. 

In Figure 1 we show with the solid line the calculated T c /T c0 as a function of 
7 7RT /T c0 (note that for 7, rcr /T c o > 2 the scale is logarithmic). The dotted line represents 
the prediction based on the BCS weak coupling approach of Ref. [19] where T c drops 
initially with the slope — n/8 (see Fig. 1 and Eq. (13) in [19]). Clearly the full strong 
coupling calculation with realistic electron-phonon spectral functions and Coulomb 
repulsion parameters leads to a much slower drop in T c with increasing interband 
impurity scattering rate than what was obtained in [19] using the BCS approach: for 
'-f 7Ta . = T c0 we get a drop in transition temperature of only 8 %, while the BCS treatment 
predicts a drop of about 36 %. 

Before we address in the next section the reasons for such a large discrepancy 
between the strong coupling and the BCS results we give in Figs. (2) and (3) calculated 
functional derivatives ST c /Sa 2 Fij, i,j = <j,tv which show how the sensitivity of 
T c to various electron-phonon couplings changes with increasing interband impurity 
scattering. In [25] these functional derivatives were computed for the case of no impurity 
scattering and it was found that the band-diagonal functional derivatives ST c /Sa 2 F aa 
and ST C / Sa 2 'F n7r have similar shapes but their overall magnitudes differ due to difference 
in sizes of the gap- functions A cr (n) and A 7r (n) in the two bands. This is also the 
case for the lowest 7^ in Fig. 2, but as the scattering rate increases the difference in 
magnitudes of the two gaps becomes smaller and the overall magnitudes of the two band- 
diagonal functional derivatives become comparable. For = 100T c0 the magnitude of 
ST C / Sa 2 F n7T is larger than the magnitude of ST C / 5a 2 F aa presumably because A W7r < \ aa 
[28] and the difference in sizes of the gaps A CT (n) and A 7r (n) has largely disappeared. 
Another consequence of this disappearance of the difference between the gaps in the 
two bands is that with increasing the shapes of 5T c /8a 2 F na and ST c /Sa 2 F aw , Fig. 3, 
become more and more similar to the shapes of the band-diagonal functional derivatives. 
The divergencies at Q = still persist, but are progressively confined to smaller and 
smaller neighborhoods of Q — 0. Again, the local maximum in ST c /8a 2 F wa near Q 
= 8k B T c for 7^0- = 100T c0 is higher than the corresponding maximum in 5T c /5a 2 F a7r 
because \ na < \ a7r [28]. 




Figure 1. The relative change in the transition temperature T c /T c0 as a function of 
the intcrband impurity scattering rate parameter "f wa -/T C Q. The solid curve gives the 
results of the full strong coupling calculation for MgB2 and the long dashed curve gives 
the results obtained from 9-9 model with interaction parameters that were used in the 
strong coupling calculation. The dotted line indicates the weak coupling result from 
Ref. [19] which predicts for low values of 7^ a straight line with the slope — tt/8. 



3. 9-9 model calculation 

The main difference between the strong coupling Eliashberg theory and the BCS 
approach is that the latter does not include the renormalization and the damping effects 
associated with the electron-phonon interaction. In the BCS calculation Zf(n), Eq. (4), 
is set equal to 1. It is possible to improve upon the BCS approach so that the effects 
of renormalization by electron-phonon interaction are included in an approximate way 
through so-called 9-9 model [26]. In this model Xij(n — m) in Eq. (2) for the electron- 
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Figure 2. The band-diagonal functional derivatives 5T c /5a 2 F aa (solid lines) and 
ST c /Sa 2 F n7r (long dashed lines) for four different values of 7 wcr given in units of T c0 . 




Figure 3. The band-off-diagonal functional derivatives 8T C / Sc^F^y (solid lines) and 
5T C / 5a 2 F a n (long dashed lines) for four different values of 7^ given in units of T c0 . 
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phonon contribution to the pairing self-energy is replaced by Xij9(Q m — \u n \)9(Q m — \u m \) 
with Q m - the maximum phonon energy (BCS approximation) and, after rewriting the 
sum in (4) as 



UJ r 



n-1 

2£ A 

m=l 



m 



Ay(0) , 



Xij(m) in Eq. (4) for the electron-phonon contribution to the renormalization function 
is replaced by Xij8(Q m — \ou m \). After rescaling the Coulomb repulsion matrix from 
cutoff uo c to Vl m according to the general prescription (11), i.e. fi*{Vl m ) — [1 + 
fi*(uj c ) l n (^ c /^m)]~ 1 /i*(^c) one gets instead of (7) 



A,(n) 



E A y + E 



1 1 



for \u n \ < Q m . 



j^i \LU n \ 

After defining 



E (Ay 

j,|w m |<n m 



,,*(0 1WT A J^ 



'i + E A ^> 



A y = 



VI + Sfc AjfcVl + Efc 



(17) 

(18) 

(19) 
(20) 



where tn = r i2 = r 2 iN 1 /N 2 , t 12 = t 21 = t 21 ^Ni/N 2 , t 22 = t 21 (in this section we label 
the a band with 1 and n band with 2), Eq. (16) can be written as 



\2n 



1 
1 



+ 



-G 2 i 



-G\ 2 
G 22 



Sin 
S 2n 



An A12 
A21 A22 



E 

|w m |<n m 



Sim 

s 2m 



(21) 



or, realizing that the right-hand side of (20) does not depend on the Matsubara index 
and denoting the elements of the corresponding 2x1 matrix with c\ and c 2 , as 

-1 / \ / \ 



An A12 
A21 A22 



E 

\u! m \<Cl r , 

The 2x2 matrix 



12771 - 11 



1 
1 



Gu 
-G 2 i 



-Gi 2 
G 22 



Cl 

c 2 



Cl 

c 2 



.(22) 



G 



(23) 



Gu —Gi 2 \ 
\ —G 2 \ G 22 J 

is a real symmetric matrix, Eq. (19), with eigenvalues d = Gu + G 22 and (see (19) 
and subsequent definitions of various tij parameters) and could be diagonalized through 
an orthogonal transformation 

^ 0\ 
,00' 



RGR- 1 



(24) 
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where the elements of R are R n = \jGn/(G n + G22), R12 = -^22/(^11 + G22), 
R21 = —R12 and R22 = Rn- Expressing G in (21) in terms of the right-hand side of 
Eq. (23) and using 

1 , / a 



E 



\u>m\<Vn 



\2m -l\ + d 



4 




where ip is the digamma function [26], Eq. (21) can be rewritten as 




(25) 



(26) 



(27) 



The transition temperature is the highest T c for which the larger eigenvalue of M is 
equal to 1. 

We have solved Eqs. (25-26) for T c as a function of the interband impurity scattering 
rate and our results are shown by the long dashed line in Fig. I. The electron-phonon 
interaction parameters were taken to be the same as those used in section 2: An = \ aa 
= 1.017, A12 = A CT7r = 0.212, A22 = A^ = 0.446 and A21 = A^ = 0.155. The maximum 
phonon energy was taken to be Q m = 75 meV, which is roughly the position of the 
largest peak in a 2 F aa (see Fig. 1 in Ref. [25]) and the values of fj,*(u c )s from section 
2.2 were scaled down using (11) to the new cutoff Q m to give //^(fim) = /x* CT (f2 m ) = 
0.139578, fi* 22 (tt m ) = /i^(O m ) = 0.139217, fJ,* 12 {Q m ) = fi* 21 (tt rn ) = 0.027081. Clearly, 
including the electron-phonon renormalization effects improves the BCS treatment 
considerably. However, we want to stress that 9-9 model gives the improved values 
only for the reduced quantity T c /T c0 as a function of the reduced interband scattering 
rate 7^0-/^0- The absolute values of T c are not accurately predicted by 9-9 model (e.g. 
we get too large a value for T c0 of 143 K so that the usual weak coupling approximation 
ip(Q m /(27iT c0 ) + l)—ip(l/2) ps ln(2e 7 f2 m /(7rT c o)), where 7 is the Euler's constant, cannot 
be made). 



4. Summary 

We have calculated the change in the superconducing transition temperature of MgB 2 
caused by interband nonmagnetic impurity scattering using the Eliashberg theory with 
realistic electron-phonon [16] and Coulomb repulsion [24] parameters for this compound. 
Our central result is given in Fig. 1. We find much slower rate of T c suppression than 
what is obtained from the BCS approach [19] which ignores the renormalization and 
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damping effects associated with the electron-phonon interactions. For small interband 
scattering rates the strong coupling calculation gives about 4.5 times slower suppression 
rate of T c than the BCS approach. Moreover, the strong coupling calculation indicates 
that it is unrealistic to expect the transition temperature of MgB 2 to ever drop below 
60% of its maximum value as a result of impurity scattering and the 20% drop in T c 
upon irradiation by fast neutrons [13] is certainly within our calculated range (20% 
drop in T c is obtained with 7 7r<T of about 4/c B T c0 , Fig. 1). Hence, the initial expectations 
based on the BCS treatment of the two-band model [19] that a dramatic suppression of 
T c in MgB 2 with interband impurity scattering would provide the final "smoking gun" 
evidence for the two-band model was exaggerated. 

Our calculation with 9-9 model (long dashed line in Fig. 1) clearly indicates that 
the main reason for the failure of the BCS approach to quantitatively account for the 
dependence of T c /T c q on r ) nc /T C Q in MgB2 is that the BCS treatment leaves out the 
electron-phonon renormalization effects. One should keep in mind, however, that for 
other multiband systems with electron-phonon and Coulomb interactions different from 
those calculated [16,24] for MgB 2 one would have to recalculate T c /T c0 as a function 
of interband scattering rate ^a/Tco using Eliashberg equations from section 2.1 with 
the interaction parameters which are relevant to the multiband superconductor that is 
being considered. 
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